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Abstract 

PaC0 2 affects cerebral blood flow (CBF) and its regulatory mechanisms, but 
the interaction between neurovascular coupling (NVC), cerebral autoregula- 
tion (CA), and cerebrovascular reactivity to C0 2 (CVR), in response to hyper- 
capnia, is not known. Recordings of cerebral blood flow velocity (CBFv), 
blood pressure (BP), heart rate, and end-tidal C0 2 (EtC0 2 ) were performed 
in 18 subjects during normocapnia and 5% C0 2 inhalation while performing 
a passive motor paradigm. Together with BP and EtC0 2 , a gate signal to 
represent the effect of stimulation was used as input to a multivariate autore- 
gressive-moving average model to calculate their separate effects on CBFv. 
Hypercapnia led to a depression of dynamic CA at rest and during stimula- 
tion in both hemispheres (P < 0.02) as well as impairment of the NVC 
response, particularly in the ipsilateral hemisphere (P < 0.01). Neither hyper- 
capnia nor the passive motor stimulation influenced CVR. Dynamic CA was 
not influenced by the motor paradigm during normocapnia. The CBFv step 
responses to each individual input (BP, EtC0 2 , stimulation) allowed identifi- 
cation of the influences of hypercapnia and neuromotor stimulation on CA, 
CVR, and NVC, which have not been previously described, and also con- 
firmed the depressing effects of hypercapnia on CA and NVC. The stability of 
CVR during these maneuvers and the lack of influence of stimulation on 
dynamic CA are novel findings which deserve further investigation. Dynamic 
multivariate modeling can identify the complex interplay between different 
CBF regulatory mechanisms and should be recommended for studies involv- 
ing similar interactions, such as the effects of exercise or posture on cerebral 
hemodynamics. 



Introduction 

Neural activation, induced by cognitive or sensorimotor 
paradigms, leads to increases in cerebral blood flow 
(CBF), a mechanism known as neurovascular coupling 
(NVC) (Girouard and Iadecola 2006). Quantification of 
the NVC response to stimulation is commonly performed 
by assessing the percent change from baseline values in 
CBF or CBF velocity (CBFv), recorded with transcranial 
Doppler ultrasound (TCD) or functional MRI, without 
taking into account the potential influence of other 



noncerebral hemodynamic variables. More recently, the 
limitations of this approach have become clear with the 
demonstration that the CBF/CBFv response is significantly 
influenced by concomitant changes in arterial blood pres- 
sure (BP) and PaC0 2 (Moody et al. 2005; Panerai et al. 
2012b; Salinet et al. 2012a). 

These observations have considerable implications for a 
better understanding of physiological and clinical studies 
of NVC. Changes in BP induced by the maneuver will 
trigger an autoregulatory response, contributing to further 
transient changes in CBFv. As shown by Aaslid and 
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coworkers (Aaslid et al. 1989; Tiecks et al. 1995), rapid 
changes in BP are followed by similar changes in CBFv, 
but due to the mechanism of dynamic cerebral autoregula- 
tion (CA), CBFv will return to its baseline value within a 
few seconds, even if BP remains displaced from baseline 
values. Moreover, changes in PaC0 2 , often estimated 
from continuous recordings of end-tidal C0 2 (EtC0 2 ), 
will also lead to changes in CBF due to the well-known 
mechanism of cerebrovascular reactivity to C0 2 (CVR) 
(Poulin et al. 1996; Brian 1998; Vernieri et al. 1999; Bat- 
tisti-Charbonney et al. 2011; Willie et al. 2012). 

In a previous study, we have shown that hypercapnia 
impairs the CBFv response to passive motor stimulation 
(Maggio et al. 2013). Analysis of instantaneous pressure- 
velocity curves indicated that the main effects resulted 
from a reduced critical closing pressure (CrCP) response, 
thus suggesting an impairment of metabolic control 
(Panerai et al. 2005, 2012a; Maggio et al. 2013) On the 
other hand, resistance-area product (RAP), which could 
be an indicator of myogenic control (Panerai et al. 2005), 
was not affected by hypercapnia, thus raising doubts 
about the influence of hypercapnia on dynamic CA dur- 
ing neural stimulation. To compare the efficacy of 
dynamic CA at rest and during neural stimulation, at 
both normocapnia and hypercapnia, using conventional 
approaches, such as transfer function analysis (Zhang 
et al. 1998; Panerai et al. 1999), is technically complex 
due to the nonstationarity induced by motor stimulation 
and also for neglecting the dynamic contribution of fluc- 
tuations in PaC0 2 . To overcome this difficulty, multivari- 
ate models have been proposed using as input variables 
both fluctuations in BP and EtC0 2 (Panerai et al. 2000; 
Edwards et al. 2004; Mitsis et al. 2006). Using a similar 
approach, we have reanalyzed the same set of data in 
(Maggio et al. 2013) using a more elaborate model 
recently described (Panerai et al. 2012b) which can pro- 
vide simultaneous estimates of dynamic CA and CVR at 
rest, together with estimates of NVC during motor stimu- 
lation. Our objective was to test the hypothesis that all 
three CBF regulatory mechanisms (CA, CVR, NVC) are 
impaired by mild hypercapnia. 

Materials and Methods 

Subjects and measurements 

This study was based in the same set of data collected by 
Maggio et al. (2013). Healthy right-handed middle-aged 
volunteers (>45 years old) were recruited from depart- 
mental staff and their relatives. Right-handedness was 
assessed with the Edinburgh inventory (Oldfield 1971). 
Exclusion criteria included physical disease in the upper 
limb, poor insonation of both temporal bone windows, 



and any history of cardiovascular, neurological, or respi- 
ratory disease. The Nottingham Research Ethics Commit- 
tee 1, United Kingdom (Ref: ll/EM/0016), approved the 
study and informed written consent was obtained from 
all participants. 

Volunteers avoided caffeine, alcohol, and nicotine for 
>12 h before attending a quiet laboratory with constant 
ambient temperature of 24°C. Beat-to-beat BP was 
recorded continuously using the Finometer® device 
(FMS, Finapres Measurement Systems, Arnhem, Nether- 
lands), attached to the middle finger of the left hand. The 
servo-correcting mechanism of the Finometer" was 
switched on and then off prior to measurements. Heart 
rate (HR) interval was recorded using a 3-lead electrocar- 
diogram (ECG) and EtC0 2 was measured via nasal 
prongs (Salter Labs) by a capnograph (Capnocheck Plus). 
Bilateral insonation of the middle cerebral arteries 
(MCAs) was performed using TCD (Viasys Companion 
III; Viasys Healthcare) with a 2 MHz probe, which was 
secured in place using a head-frame. Hypercapnia was 
induced by the inhalation of a mixture of 5% C0 2 /air 
through a mask. During C0 2 breathing, the capnograph 
was connected to the mask via a sample line. During the 
entire procedure, subjects were in a supine position and 
detailed instructions were given before taking measure- 
ments. 

After a period of 15-min stabilization, participants per- 
formed a 5-min baseline recording and two passive motor 
paradigms during air breathing. The same sequence was 
repeated during C0 2 inhalation. Arterial BP was mea- 
sured with a sphygmomanometer before each measure- 
ment. 

The paradigm was performed only with the dominant 
arm and consisted of an examiner performing repetitive 
flexion and extension of the subject's elbow within a 
range of movement of approximately 90° at a rate of 
1 Hz, given by the sound of a metronome. Subjects were 
instructed to relax and not move the arm. All paradigm 
recordings started with a 90 sec baseline phase. Thereaf- 
ter, the paradigm was performed over 60 sec, with a 
90 sec recovery phase. During the rest and recovery peri- 
ods, the examiner kept hold of the participant's arm. 

Data analysis 

Data were simultaneously recorded onto a data acquisi- 
tion system (PHYSIDAS, Department of Medical Physics, 
University Hospitals of Leicester) for subsequent off-line 
analysis. ECG, EtC0 2 , BP, and stimulus marker signals 
were sampled at 500 samples s - , and BP was calibrated 
at the start of each recording. All signals were visually 
inspected to identify artefacts and noise, and narrow 
spikes (<100 ms) were removed by linear interpolation. 
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The CBFv channels were subjected to a median filter and 
all signals were low-pass filtered with a cut-off frequency 
of 20 Hz. The R-R interval was then automatically 
marked from the ECG and continuous HR plotted against 
time. Occasional ectopic beats caused spikes in the HR 
signal; these were manually removed by remarking the 
R-R intervals for the time points at which they occurred. 
Mean BP and CBFv values were calculated for each car- 
diac cycle. The end of each expiratory phase was detected 
in the EtC0 2 signal, linearly interpolated, and resampled 
with each cardiac cycle. Beat-to-beat data were spline 
interpolated and resampled at 5 samples s~ to produce 
signals with a uniform time-base. With the use of the 
electrical output from the metronome, the stimulation 
signal s(t) was added to the ensemble. The maneuver that 
achieved the highest amplitude of contralateral CBFv 
response was chosen to represent the participant's 
response (Salinet et al. 2012b). 

A multivariate autoregressive-moving average (ARMA) 
model was used to represent the dynamic influence of 
BP, EtC0 2 , and stimulus [s(t)] (inputs) on CBFv (out- 
put). These time-domain models have been used previ- 
ously to analyze the simultaneous effects of BP and 
PaC0 2 on CBFv (Panerai et al. 2000; Edwards et al. 2004; 
Mitsis et al. 2004, 2006). More recently, we have extended 
their application to NVC by adding the effects of stimula- 
tion through the s(t) input (Peebles et al. 2012; Salinet 
et al. 2012b). The order of these models, representing the 
number of past samples adopted for the autoregressive 
(AR) and moving average (MA) terms, was thoroughly 
considered (Peebles et al. 2012). As described in detail 
elsewhere (Peebles et al. 2012; Salinet et al. 2012b), 
ARMA coefficients were used to estimate the CBFv step 
response to each individual input which then represents 
the efficacy of the three separate regulatory mechanisms, 
that is CA, CVR, and NVC. Using the model autoregres- 
sive and moving average coefficients, it is possible to 
obtain the CBFv step response to each of the separate 
inputs. Step responses for each subject were averaged to 
obtain the population mean and SD step responses. For 
the CBFv step response for the BP input, which represents 
dynamic CA (Panerai et al. 2012b), the autoregulation 
index (ARI) proposed by Tiecks et al. (1995) was 
extracted by using the best least-squares fit between the 
CBFv step response and one of the 10 model ARI curves 
proposed by Tiecks et al. (1995). ARI was computed for 
each subject separately for left and right hemisphere dur- 
ing normocapnia and hypercapnia. The CBFv step 
responses to the EtC0 2 and s(t) inputs were quantified by 
their change from baseline to the plateau region, defined 
as the mean value of the response from 20-25 sec, 
expressed by parameters AS C o2 an d AS s(f ), respectively. 
AS C02 is a measure of CVR to C0 2 . During baseline, the 



same multivariate model was adopted to estimate the 
CBFv step responses to the BP and EtC0 2 inputs, by 
setting the third input (s(t)) to zero and corresponding 
values of ARI and AS C o2 were computed as for the stim- 
ulation phase. 

Statistical analysis 

Mean values of each variable were extracted from the 30 sec 
phase before the start of the paradigm to reflect baseline 
conditions. Student's T-test for dependent variables was 
used to compare baseline values of CBFv, BP, HR, EtC0 2 , 
and ARI between normocapnia and hypercapnia. 

To compare changes in CBFv between normocapnia and 
hypercapnia, the area-under-the-curve (AUC) was calcu- 
lated for their differences from the beginning of the maneu- 
ver, up to 20 sec after the end of passive arm movement. 
Statistical analysis was performed using two-way repeated 
measures ANOVA with type (normocapnia vs. hypercap- 
nia) and side (contralateral vs. ipsilateral) as within factors 
for CBFv variations. Two-way repeated measures ANOVA 
was also adopted to test for differences in ARI and AS C o2 
due to hypercapnia and the effect of the maneuver in com- 
parison to baseline values. When appropriate, a post hoc 
test for multiple comparisons (Tukey's test) was performed. 
A value of P < 0.05 was adopted to indicate statistical sig- 
nificance. All statistical analyses were performed with Statis- 
tica® software for Windows. 

Results 

Nineteen participants (eight male) of mean age 58.1 years 
(SD 7.1, range 48-77) were recruited. Participant's mean 
Edinburgh Inventory for right-handedness was 96.5% 
(8.5). All participants completed normocapnia measure- 
ments. Three participants performed the baseline record- 
ing and only one passive motor maneuver with 5% C0 2 
due to intolerance to the mask. In these cases, we ana- 
lyzed the only available passive motor paradigm recording 
to represent the participant's response. 

As described previously (Maggio et al. 2013), hypercap- 
nia increased CBFv from 56.6 (11.8) to 72.0 (14.8) 
cm s -1 (P < 0.001), in the contralateral MCA, and from 
49.4 (9.1) to 62.2 (11.2) cm s _1 (P < 0.001) in the ipsilat- 
eral MCA. Mean BP increased from 99.4 (11.1) to 109.1 
(19.3) mmHg (P < 0.01) and EtC0 2 increased from 40.7 
(3.3) to 47.9 (3.2) mmHg (P < 0.001). 

CBFv responses and its contributors 

One subject was excluded from further analysis, as he did 
not show a CBFv response to the paradigm during 
normocapnia. 
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In both normocapnic and hypercapnic conditions, 
CBFv increased bilaterally with the beginning of passive 
movement and showed a plateau phase that outlasted the 
duration of stimulation (Fig. 1). This characteristic bilat- 
eral rise in CBFv (AUC 105.2 (4.6) % ipsilateral; 104.9 
(3.8) % contralateral) was significantly reduced by hyper- 
capnia (AUC 102.2 (3.0) % ipsilateral; 103.2 (2.9) % con- 
tralateral, ANOVA, P = 0.01). A post hoc comparison 
indicated a significant difference for the ipsilateral hemi- 
sphere only (Tukey, P = 0.001). 

Bilateral CBFv step responses to the BP input extracted 
during baseline and motor stimulation are represented in 
Figs. 2A,B and 3A,B, respectively. From these step 
responses, ARI was calculated during normocapnia and 
hypercapnia. At baseline, ARI values decreased during 
hypercapnia from 6.1 (1.0) to 4.8 (1.1) (P < 0.002) in the 
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Figure 1. (A) Representative CBFv response from a single subject 
during normocapnia (continuous line) and hypercapnia (dashed line) 
to repetitive elbow flexion from 0- 60 sec (gray bar). Population 
averages of CBFv responses for the ipsilateral (B) and contralateral 
(C) middle cerebral arteries. Error bars represent the largest ±1 SE. 



ipsilateral hemisphere, and from 6.1(1.3) to 4.8 (1.9) 
(P < 0.01) in the contralateral hemisphere. Two-way 
ANOVA revealed significant differences between groups 
(F = 4.1, P = 0.001), and the post hoc Tukey test showed 
a significant decrease in ARI values during hypercapnia 
(P = 0.002 and 0.01 for the contralateral and ipsilateral 
hemispheres, respectively). During stimulation, the mean 
ARI values decreased from 6.1 (0.7) in the contralateral 
hemisphere and 6.1 (1.2) in the ipsilateral during normo- 
capnia to 5.0 (1.6) and 5.2 (1.2) during hypercapnia, 
indicating a significant impairment in CA (Tukey 
P = 0.005 and P = 0.02 for the contralateral and ipsilat- 
eral hemispheres, respectively). The similarity between 
these values and the corresponding ARI during baseline 
was confirmed by the result of the two-way ANOVA 
showing that the passive motor maneuver (during nor- 
mocapnia) had no significant effect on dynamic CA. 

Corresponding step responses for the EtC0 2 (during 
baseline and stimulation) and motor stimulus inputs 
showed similar patterns, rising gradually to reach a pla- 
teau within 15 sec (Figs. 2C,D and 3C,D). The CBFv step 
response to C0 2 and the ASco2 parameter were not sig- 
nificantly changed by either hypercapnia or passive elbow 
flexion. At baseline, the mean values of ASco2 were 2.1 
(1.3) cm s~' mmHg- 1 and 1.5 (1.6) cm s 1 mmHg 1 for 
the ipsilateral hemisphere, and 1.4 (1.8) cm s~' mmHg -1 
and 0.8 (1.7) cm s~' mmHg -1 for the contralateral hemi- 
sphere during normocapnia and hypercapnia, respectively. 
During stimulation, ASco2 was 2.8 (1.6) cm s _1 mmHg -1 
and 2.6 (0.8) cm s _1 mmHg -1 ipsilateral and 1.6 (2.0) 
cm s 1 mmHg- 1 and 1.4 (1.9) cm s 1 mmHg 1 contra- 
lateral, again for normocapnia and hypercapnia, respec- 
tively. For the CBFv step responses to the stimulus input 
(Fig. 3E,F), AS s(t , varied from 0.3 (0.09) (A.U.) and 0.4 
(0.08) during normocapnia to 0.09 (0.03) (P = 0.001) 
and 0.1 (1.0) (P = 0.04) during hypercapnia for the ipsi- 
lateral and contralateral hemispheres, respectively. 

Discussion 
Main findings 

To our knowledge, this is the first study to model the 
simultaneous influences of hypercapnia on NVC, dynamic 
CA, and CVR. Our results confirmed that dynamic CA 
was impaired by hypercapnia both at rest and during 
neural stimulation, and that NVC was also impaired by 
breathing 5% C0 2 in air. On the other hand, CVR was 
not impaired by hypercapnia or neural stimulation and 
dynamic CA was not affected by passive elbow flexion 
during normocapnia. The CBFv step responses to each 
input (BP, EtC0 2 , s(t)) had temporal patterns compatible 
with previous studies. For the BP input, the CBFv step 
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Figure 2. Population average CBFv step responses (SR) due to BP (A, B) and EtC0 2 (C, D) at baseline. Time f = 0 represent the beginning to 
the CBFv SR to each of the two inputs. Separate responses given for the ipsilateral (A, C) and contralateral (B, D) MCA during normocapnia 
(continuous line) and hypercapnia (dashed line). Error bars represent the largest ±1 SE. 



response showed the classical rapid rise, followed by a 
slower decline as expected of a working dynamic CA 
(Tiecks et al. 1995; Panerai et al. 2001). As expected, 
hypercapnia shifted the response upwards (Panerai et al. 
1999) and reduced the value of the ARI parameter (Kat- 
sogridakis et al. 2013). Of considerable interest, this 
response was not influenced by passive motor stimulation 
itself. 

The CBFv step response to the EtC0 2 input had a tem- 
poral pattern showing a gradual rise as obtained in direct 
measurements of step changes in PaC0 2 (Poulin et al. 
1996) as well as in previous studies based on system iden- 
tification techniques (Edwards et al. 2004; Mitsis et al. 
2004; Salinet et al. 2012a). The ASco2 parameter is equiv- 
alent to the classical assessment of CVR using the differ- 
ence in CBFv from baseline to a steady-state value after 
breathing a C0 2 -rich mixture with air. This parameter 
was in very good agreement with corresponding values of 
CVR published in the literature (Ide et al. 2003; Moody 
et al. 2005; Willie et al. 2012), and it was not significantly 
influenced by hypercapnia during baseline or by passive 
motor stimulation. The third step response, that quanti- 
fies the influence of stimulation on CBFv, also had a 
gradual rise to a plateau which was significantly reduced 
by hypercapnia, as shown by the complete step response 
curves (Fig. 3), as well as by the AS s(t ) parameter. These 
curves are in excellent agreement with results obtained 
from a different group of healthy subjects (Panerai et al. 



2012b), but their sensitivity to PaC0 2 had not been 
reported previously. 

Complexity of the NVC response to elbow 
flexion 

One important limitation of NVC studies based on imag- 
ing techniques is the difficulty of recording simultaneous 
changes in BP and other variables that are also disturbed 
by cognitive or sensorimotor paradigms. As a result, most 
imaging studies tend to ascribe all changes in CBF to the 
stimulus when in fact there is substantial contribution 
from peripheral determinants as shown by previous stud- 
ies (Moody et al. 2005; Panerai et al. 2005, 2012b; Salinet 
et al. 2012a,b, 2013a,b; Maggio et al. 2013; Nogueira et al. 
2013). These results cannot be generalized though as dif- 
ferent paradigms are likely to have variable degrees of 
influence on BP and autonomic nervous system control, 
as well as on respiratory frequency and other physiologi- 
cal covariates (Moody et al. 2005; Salinet et al. 2012a). 
For this reason, studies of NVC should ideally include 
simultaneous measurements of BP, EtC0 2 , and other 
physiological variables that could influence the CBFv 
response under different circumstances. Due to parallel 
changes in BP and PaC0 2 , a clear interaction exists 
between dynamic CA, CVR, and NVC. This raises the 
question of whether multivariate dynamic modeling 
should be the tool of choice in the context of NVC 
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studies, to allow for simultaneous assessment of dynamic 
CA and CVR. 

Changes in BP during C0 2 stimuli have direct effects 
on CBFv (Claassen et al. 2007), so that only in the CBFv 
response range below the threshold for the increase in BP 
with C0 2 the CBFv measurement reflect vascular reactiv- 
ity to C0 2 alone (Battisti-Charbonney et al. 2011). Simul- 
taneous modeling of the influences of BP and PaC0 2 on 
CBFv have been described previously, mainly for record- 
ings performed at rest (Peebles et al. 2012) but also for 
perturbations in PaC0 2 (Edwards et al. 2004) and ortho- 
static challenges induced with lower body negative pres- 
sure (Mitsis et al. 2006). In particular, Edwards et al. 
(2004) have demonstrated the sensitivity of CBFv step 
responses to BP and EtC0 2 for both hypo- and hypercap- 
nia, and Mitsis et al. (2004) have identified a dynamic 
interaction between the BP and EtC0 2 inputs. Neverthe- 
less, none of these models have addressed the interaction 
between dynamic CA, CVR, and NVC mechanisms. One 



recent study though has introduced a technique to com- 
pensate for the influence of BP on NVC responses 
(Spronck et al. 2012), based on Tiecks et al. (1995) step 
response templates for the ARI parameter. 

One important aspect of the interaction of different 
CBF regulatory mechanisms is whether CA and/or CVR 
are depressed during neural activation. Applying classical 
techniques of assessment, dynamic CA was found to be 
impaired during cognitive stimulation based on both 
word construction and a 2D puzzle, which preferentially 
activate the left and right hemispheres, respectively (Pane- 
rai et al. 2005). Similar conclusions could be drawn from 
other studies (Nakagawa et al. 2009; Dineen et al. 2010; 
Nogueira et al. 2013). The finding that the CA response 
was not altered by passive motor stimulation in the pres- 
ent study could be explained by differences in the inten- 
sity of the alert reaction induced by other paradigms. An 
alternative explanation though would be the fact that 
multivariate modeling takes into account the influences of 
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PaC0 2 and the stimulus, thus rendering more accurate 
estimates of dynamic CA. Similar concerns would apply 
to studies of cerebral hemodynamics during exercise, 
where CA was found to be normal during mild exercise 
but depressed during exhaustive exercise (Ogoh et al. 
2005). 

Unaccounted for interactions with BP and stimulation 
could also explain why the step responses to C0 2 (Fig. 3) 
and the AS C0 2 parameter did not show changes with 
hypercapnia per se or due to the s(t) input. The CBFv- 
arterial C0 2 relationship has been described as exponen- 
tial by some authors (Markwalder et al. 1984; Ide et al. 
2003), but others have reported a more linear relationship 
(Claassen et al. 2007; Willie et al. 2012), mainly in the 
range of 40-50 mmHg as was the case of our study. The 
latter would explain a constant CVR with increasing levels 
of PaC0 2 as we have found. Moreover, the exponential 
relationship was established with steady-state changes in 
PaC0 2 , and did not take into account interactions with 
BP (Dumville et al. 1998). However, using the two-breath 
technique, Edwards et al. (2004) obtained step responses 
to EtC0 2 changes that were depressed by hypercapnia. 
Those investigators took into account influences of BP by 
using cerebrovascular conductance, instead of CBFv, as 
the model output variable. The lack of a significant effect 
of hypercapnia on the step responses to C0 2 in our case 
could be due to the much smaller signal-to-noise ratio of 
the EtC0 2 input, resulting from spontaneous fluctuations, 
as compared to the relatively larger changes induced by 
the two-breath technique with inhalation of 8% C0 2 . 

The interaction of CVR with passive motor stimulation 
also deserves attention. The influence of sympathetic 
stimulation on classical measures of CVR is still contro- 
versial with studies showing both significant (Peebles 
et al. 2012) and nonsignificant influences (LeMarbre et al. 
2003). An increase in CVR was reported during whole- 
body exercise (Ogoh et al. 2005) and it is possible that 
the absence of an effect in our case was due to the much 
smaller muscle mass involved in passive elbow flexion. 

Overall, the considerations above suggest that in the 
presence of physiological changes induced by cognitive/ 
sensorimotor paradigms, posture, exercise, breathing, or 
other maneuvers, multivariate modeling should be the 
tool of choice. Nevertheless, further work in this area is 
needed to optimize the use of multivariate modeling in 
studies of cerebral hemodynamics. 

Influence of C0 2 on NVC 

The marked reduction in the NVC response induced by 
hypercapnia, as reflected by the CBFv step responses to 
the s(t) input (Fig. 3) and the AS s(f ) parameter, is in 
agreement with previous studies in humans and animals 



(Xu et al. 2011; Kennerley et al. 2012), but at odds with 
similar studies based on visual stimulation and CBFv 
measurements in the posterior cerebral artery (PCA). 
Rosengarten et al. (2003) found no differences due to 
hypercapnia in the CBFv response to a reading para- 
digm. On the other hand, using the same protocol, 
Szabo et al. (2011) reported a reduction in the CBFv 
response to visual stimulation in healthy young subjects 
during hypocapnia induced by hyperventilation. These 
apparently contradictory results might be explained by 
regional differences in NVC between the visual and 
motor cortex (Chiarelli et al. 2007). Although regional 
differences in CVR have been reported, based on MRI 
measurements (Ito et al. 2000), no differences in CVR 
were obtained comparing transcranial Doppler record- 
ings of the MCA with those from the PCA (Willie et al. 
2012). One important alternative possibility though is 
that the effects of changes in PaC0 2 on NVC, and dif- 
ferences between motor and visual areas, are mainly 
related to the direct influence of C0 2 on the neural 
response to stimulation, instead of the vasomotor effects 
resulting from changes in pH (Yoon et al. 2012). This 
interpretation could also explain the 'flow paradox' 
reflected by the baseline changes in CBFv when com- 
pared to the amplitude of the NVC response. In our 
study, hypercapnia increased baseline CBFv by more 
than 25% in both hemispheres. On the other hand, neu- 
ral stimulation resulting from elbow flexion, led to an 
increase in CBFv of the order of 5% (Fig. 1). Therefore, 
during hypercapnia there would be enough excess blood 
flow to satisfy the metabolic demand of neural stimula- 
tion. Nevertheless, a response could still be detected, 
albeit of lesser amplitude (Fig. 1). These observations 
suggest that the flow response is directly linked to neu- 
ral signaling and not to the feedback effect of increased 
demand for 0 2 and that it is the forward neural 
response that is affected by hypercapnia (Yoon et al. 
2012). Further support for this interpretation can be 
derived from the study of Xu et al. (2011) in humans 
showing that hypercapnia has a suppressive effect on 
brain activity, including a reduction in functional con- 
nectivity. In rats, Kennerley et al. (2012) also demon- 
strated that the differences observed due to hypercapnia 
on the response to whisker stimulation were mainly 
reflecting changes in neural activity than in the demand 
for 0 2 . The lack of effect of hypercapnia during visual 
stimulation (Rosengarten et al. 2003) suggests that the 
direct influence of C0 2 on neural activity might be dif- 
ferent in the visual cortex and could also explain the 
results of hypocapnia induced by hyperventilation (Szabo 
et al. 2011). More work in this area should be a prior- 
ity, including studies of the effects of hyperventilation 
during sensorimotor paradigms. 
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Limitations of the study 

Changes in CBFv can only reflect corresponding changes 
in CBF if the diameter of the insonated artery remains 
constant, an assumption that has been supported by 
observational and imaging studies (Giller et al. 1993; Ser- 
rador et al. 2000). These studies did not show changes in 
MCA diameter due to hypercapnia, but similar investiga- 
tions are lacking for the potential influence of neural 
stimulation. If MCA diameter should increase due to pas- 
sive elbow flexion, our recorded changes in CBFv due to 
stimulation would be an underestimate of the true 
changes in CBF thus biasing our results in the direction 
of type II errors. 

The multivariate model, based on ARMA structures is 
limited to a linear representation of input-output rela- 
tionships. Although nonlinear models have been proposed 
for dynamic CA, their superiority has not been demon- 
strated (Panerai et al. 2001; Mitsis et al. 2004). Neverthe- 
less, one aspect of multiple input nonlinear models is the 
possibility to incorporate input-input interactions, as 
demonstrated by Mitsis et al. (2004) for the BP and 
EtC0 2 inputs. Further work is needed to explore the rele- 
vance of incorporating this type of interaction in the 
model, mainly from the perspective of improving the reli- 
ability of clinical applications. 

Separate training and testing sets of data were not 
appropriate to the nature of this study as the main objec- 
tive was not to propose a general model for CA, CVR, 
and NVC, but simply to use the flexibility of the ARMA 
structure to extract parameters leading to the estimation 
of CBFv step responses to each input. This is equivalent 
to the application of linear regressions to extract the slope 
of static relationships. 

Motor stimulation with the dominant arm would be 
expected to lead to pronounced interhemispherical differ- 
ences. In our study, contralateral CBFv step responses to 
stimulation (Fig. 1) were significantly greater in compari- 
son with the ipsilateral side thus confirming the physio- 
logical consistency of the model. On the other hand, 
lateralization was somewhat limited, possibly due to the 
nature of the paradigm, involving more diffuse and non- 
specific stimulation regarding attention, concentration, 
and motivation. 

Finally, in our previous communication (Maggio et al. 
2013), ARI values were smaller than in the current study. 
In both cases, ARI values were estimated by fitting Tiecks 
et al. (1995) template responses, but in our first study, 
CBFv step responses were obtained by FFT-based transfer 
function analysis, while in the current study, step 
responses were estimated using ARMA models. In previ- 
ous studies of dynamic CA, we also observed that mean 
normal values obtained with ARMA modeling were 



higher than those obtained with classical transfer function 
analysis (Panerai et al. 2003, 2010). 

Conclusions 

Hypercapnia induces cerebral vasodilation and impair- 
ment of CBF regulatory mechanisms, chiefly CA. We con- 
firmed that NVC was also impaired by breathing 5% C0 2 
in air, but CVR was not influenced either by hypercapnia 
or neural stimulation, and dynamic CA was not influ- 
enced by passive motor stimulation. The NVC response 
during hypercapnia can be better explained by a direct 
effect of C0 2 on neural activity, rather than the vasodila- 
tion associated with reduced pH. More research is needed 
on the lack of change in CVR due to hypercapnia. Neural 
stimulation by cognitive or sensorimotor paradigms 
induces simultaneous changes in BP and PaC0 2 , thus 
requiring the use of multivariate dynamic models for 
identification of its separate influences and more robust 
assessment of the NVC response. Further application of 
this approach is warranted under diverse physiological 
and clinical conditions. 
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